Mon. Not. R. Astron. Soc. 000, [j-ju] (2001) Printed 1 February 2008 (MN I*TeX style file vl.4) 



Convective cores in galactic cooling flows 



Alexei Kritsuk^'^, Tomasz Plewa^'^ & Ewald Miiller^ 

^ Max-Planck-Institut fiir Astrophysik, Postfach 1317, D-85741 Garching, Germany 

^ Sobolev Astronomical Institute, University of St Petersburg, Stary Peterhof, 198904 St Petersburg, Russia 
Nicolaus Gopernicus Astronomical Center, Party cka 18, 00716 Warsaw, Poland 



Submitted: August 1, 2000; Revised: March 21, 2001 



o 
o 

(N 
>> 



m 
> 

o 

oo 
O 

o 
o 



6 



ABSTRACT 

We use hydrodynamic simulations with adaptive grid refinement to study the depen- 
dence of hot gas flows in X-ray luminous giant elliptical galaxies on the efficiency of 
heat supply to the gas. We consider a number of potential heating mechanisms includ- 
ing Type la supernovae and sporadic nuclear activity of a central super-massive black 
hole. A s a starting po int for this research we use an equilibrium hydrostatic recycling 
model ( Kritsuk 1996 ). We show that a compact cooling inflow develops, if the heating 
is slightly insufficient to counterbalance radiative cooling of the hot gas in the cen- 
tral few kiloparsecs. An excessive heating in the centre, instead, drives a convectively 
unstable outflow. We model the onset of the instability and a quasi-steady convective 
regime in the core of the galaxy in two-dimensions assuming axial symmetry. 

Provided the power of net energy supply in the core is not too high, the con- 
vection remains subsonic. The convective pattern is dominated by buoyancy driven 
large-scale mushroom-like structures. Unlike in the case of a cooling inflow, the X-ray 
surface brightness of an (on average) isentropic convective core does not display a 
sharp maximum at the centre. A hybrid model, which combines a subsonic peripheral 
cooling inflow with an inner convective core, appears to be stable. We also discuss 
observational implications of these results. 
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1 INTRODUCTION 

Extended thermal X-ray emission from the hot (~ 10^ K) in- 
terstellar medium (ISM) in giant elliptical galaxies is usually 
interp reted in terms of a 'cooling flow' scenario [see Loewen- 
stein (1997) for a recent review]. It implies that radiative 



cooling of optically thin hot plasma drives a subsonic inflow 
towards the centre of a potential well formed by the stellar 
component and a massive dark halo. Galactic cooling flows 
and cooling flows in clusters of galaxies are essentially dif- 
ferent owing to a difference in the origin of the two media. 
The intracluster medium is mostly primordial, cools from 
its high virial temperature and is accreted by the central 
galaxy, which provides a focus for the flow. In contrast, the 
ISM in elliptical galaxies is constantly replenished by mass 
loss from evolved stars. This gas must be thermalized in 
the galaxy's gravity field and heated to X-ray temperatures 
from which it may cool down again. Thus, in hydrodynamic 
terms, galactic cooling flows are deflned as systems where 
the mass flux is dominated by source terms a s opposed to 
'boun dary terms' as in cluster cooling flows (Loewenstein 
199^). 



counterbalance radiative cooling of the hot ISM in ellipticals. 
The energy supplied with normal stellar mass loss is limited 
by the value of the stellar velocity dispersion and may only 
provide tempe ratures a factor of ~ 2 lower than the ISM 
temperatures (Davis & White 1996). In contrast, gravita- 



tional energy liberated in case of gas inflow wou ld pr ovide 
too much heat to the central regions [Thomas (L986) and 
references therein]. 

Heating by supernovae (SN) Type la remains a contro- 
versial issue. On the one hand, the SN rate is not very well 
constrained. Its value varies from 0.55 /175 SNu, suggested by 



Several potential heat sources have been considered to 



van den Bergh & Tammann ( |l99l| ), the uncertainty being a 
factor of the order of 1.5, to 0. 18 ± .06 h^^ SNu derived by 
Cappellaro, Evans & Turatto ( |l999| ).p| 

On the other hand, the low iron content of the diffuse 
hot gas in ellipticals estimated from high-quality broad band 
X-ray spectra provided by ASCA for a single-t emperature 
plasma model (Loewenstein & Mushotzky 199J) suggests a 
very low efficiency of ISM enrichment by Type la SNe. This 



/i75 is the Hubble constant in units of 75 km s ^Mpc ^\ 1 SNu 
: 1 supernova per century per 10^" solar bolomctric luminosities. 
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implies either that the SN rate is lower than ~ 0.03 /175 SNu 



( Loewenstein fc Mushot zky 199^ ) or t hat SN ejecta do not 
mix wi th the hot ISM (Thomas 19St; Fujita, Fukumoto & 
Okoshi 1997). However, the multi-parametric analysis of X- 
ray spectra is a complex procedure based on iron L line diag- 



nostic t ools, and requires accurate atomic physics data ( Ari- 
moto e . al. 1997). The procedure is also model-dependent 



Using two-temperature multi-phase plasma models Buote 
( 1999| ) obtained substantially better spectral fits for the 
same data sets with iron abundances of ~ 1 — 2 solar and 
relative element abundances fixed at their solar values. His 
results are consistent with the Type la SN rate being up 
to a factor of ~ 2 lower than reported by Cappellaro et 



al. ( 199£ ). Clearly, better quality data are required to deter- 
mine the Type la SN heating rate more precisely. 

All of the above mentioned energy sources may be de- 
scribed as being continuously distributed within a galaxy 
and their local rates depending on the stellar mass density 
and velocity dispersion, the gas inflow velocity profile, and 
the shape of the gravita tional potential. There is a 'global 
problem' (Thomas 198^) of balancing heating and cooling 
both in the centre and in the outer regions, since the source 
terms depend in different ways on physical characteristics 
which vary with radius. However, empirical scaling laws for 
ellipticals and the physics of radiative cooling and thermal 
instability imply certain restrictions on the radial depen- 
dences. One possible solu tion to this pr oblem is a hydrostatic 
hot gas recycling model ( Kritsuk 1996 ) for hot coronae of el- 
liptical galaxies. In the hydrostatic gas configuration all that 
remains of the gas dynamical equations are the algebraic 
source terms which balance to zero. In this model two equi- 
librium conditions are simultaneously satisfied: (1) the stel- 
lar mass loss rate exactly balances the rate that mass cools 
locally from the fiow (dropout) and (2) the rate that thermal 
energy is radiated away is exactly balanced by the energy 
input from stars and supernovae. The recycling model can 
be used as a tool to distinguish between infiow and out- 
flow regimes for a given set of galaxy parameters (Kritsuk, 
Bohringer & Miiller 1998). 

Periods of sporadic activity can drastically change the 
thermal state of the ISM in the inner 'cooling flow' region 
of a cluster or of an elliptical on a time scale of ~ 10* yr. 
There are numerous observations of an interaction between 
radio sources and the hot gas both in centra l cluster galax- 
ies and o ther ellipticals (iHarris et al. 1999 ; [Koekemoer et 
al. 199^; |McN amara et al. 2000[ ). It is this complex central 
region that is most sensitive to any imbalance in the energy 
budget of the ISM. Since the 'thermal' time scale is shorter 
there, the core of a few kpc in radius is the first to reach 
a new equilibrium state corresponding to the current status 
of heating. High spatial and spectral resolution X-ray ob- 
servations of these central regions are critical to distinguish 
among the various heat sources and flow regimes in the core. 

The response of a cooling flow to energy injection by 
the central black hole that is fed by the cooling flow was 
simulated by Binney & Tabor (1995). In their spherically 
symmetric numerical models the black hole heats the cen- 
tral kiloparsec of the ISM as soon as it begins to swallow 
gas from the ISM. The resulting expansion of the heated gas 
eliminates the cusped density profile and establishes a core 
that is roughly 2 kpc in radius. After the central heat source 
has switched off, the core again cools catastrophically within 



0.5 Gyr. Among several limitations of the model (which in- 
clude a peculiar feedback mechanism) the authors particu- 
larly note the restriction of spherical symmetry. Both the jet 
and the ensuing convective instabilities cannot be properly 
modelled due to this symmetry restriction. 

In this paper we analyze the hot gas flow regimes, which 
develop in elliptical galaxies in response to a small imbal- 
ance in the energy budget caused b y a global instabilit y 
of the equilibrium recycling model (Kritsuk et al. 1998). 



We use multi-dimensional numerical simulations to inves- 
tigate compact cooling inflows and subsonic convection - 
two alternative saturated nonlinear regimes for a slightly 
disturbed thermal equilibrium in the core. Here, we restrict 
ourselves to spherically symmetric (on the average) convec- 
tive flows, which have not been studied before. The symme- 
try restriction will be relaxed in a subsequent paper which 
more closely addresses the astrophysical aspects of the prob- 
lem. The analysis presented here provides a basis for a better 
understanding of more complex convective flows and allows 
for a direct comparison with one-dimensional models. 

A short description of our basic model is given in Sec- 
tion 2. Section 3 describes the numerical method. The results 
of numerical experiments are presented in Section 4 together 
with the detailed properties of the inner convective cores. In 
Section 5 we summarize the results and discuss effects to be 
expected for the dynamical state of the hot gas in ellipticals 
from a larger imbalance between cooling and heating. 



2 THE MODEL 

Our numerical experiments are based on a physical model, 
which has been described in Kritsuk et al. (1998). We refer 



the reader to that paper and references therein for a more 
detailed description than given below. 

We assume that the local mass budget of the hot ISM 
inside the galaxy is determined exclusively by mass supply 
from the old stellar population of the galaxy [in the form of 
stellar winds, planetary nebulae, supernova explo sions , etc.; 
a detailed discussion c an be found in Mathews (1990) and 
Mathews & Brighenti (1999)], and by gas condensation to 
a warm phase due to thermal instabilities. Our motivation 
for including mass dropout is based on the fact that the 
assumed mass supply mechanisms inevitably cause nonlinear 
small scale perturbations in the hot phase, which can trigger 
thermal instabilities. The rate of mass dropout is assumed 
to be a function of local conditions in the hot ISM, and 
given by the expression pti ~ bx{n)p, where p is the density 
of the hot phase. The spectrum of the initial small-scale 
perturbations in the gas and, in particular, its modification 
by heat conduction are described by a single parameter b = 
const G [0, 1]. The Heaviside function x ^-nd the linear 
instability growth rate n are defined as 



n, if n > 0; 
0, otherwise; 
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Here T is the hot gas temperature, A(r) is the nonequilib- 
rium zero field radiative cooling function for solar compos 



tion and complete ionization of the plasma (Sutherland & 
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Dopita 1993), and Cp is the specific lieat at constant pres- 
sure. The first term in the rhs of equation coincides with 
Field's instability criterion, and the s econd one de scribes 
stabilization due to stellar mass loss dKritsuk 19*9^ ), Q IS 
the specific stellar mass loss rate, and p, is the local stellar 
mass density. Here we assume b « 0.5 (Kritsuk 1996). Since 



we are interested in relatively short evolutionary time scales 
(< 1 Gyr), we assume a to be a constant. 

The local energy budget of the hot ISM is determined 



by radi lative losses and energy supply from stars (via SN 
Type la explosions and via thermalization of the kinetic en- 
ergy of orbital motion for the gas shed by stars). Radiative 
cooling comprises two energy sinks due to (i) direct radia- 
tive cooling of the hot phase and (ii) due to radiative losses 
by the condensing material. 

Finally, the hydrodynamic description for the hot phase 
is based upon the assumption that the sources are homoge- 
neously distributed. Note that this assumption is fulfilled as 
long as the numerical resolution is not excessively high. The 
set of conservation laws for mass, momentum, and energy of 
the hot phase can then be written as follows: 



dp 
dt 

d{pv) 



dt 



+ V-(pt;) = «p* - pti, (3) 
-f uV-(pi;) -f pv ■ Vd -|- Vp = pV(j) — ptiV, (4) 



dE 
'dt 



+ 'V-[v{E+p)] = 



Pti{E + p)/p - p^K + pvV(t 



(5) 



Here v — {u, v, w) is the velocity vector, p ~ {'y — l)pe is the 
pressure, 7 is the adiabatic index, e is the specific internal 
energy, and e* — c^Tq is the characteristic specific energy 
of the heat source of temperature Tq. The energy density 
E = p{e + v''/2). 

We define a spherically symmetric gravitational poten- 
tial (f) to be that of a two-component isot hermal sphere , 
which includes stars and a massive dark halo (Kritsuk 1997). 
The galaxy model is characterized by four parameters: the 
stellar velocity dispersion a* , the characteristic radius of the 
stellar mass distribution 



ro = 



a* 



y^4:TvGp,,o 



(6) 



the ratio of the stellar and dark matter (DM) velocity dis- 



persions 



/9: 



< 1, 



(7) 



and the ratio of the dark matter mass density and the stellar 
mass density at the centre of the galaxy 



S = Pdm,o/p*,o- 



(8) 



For the prototype galaxy NGC 4472, we set ct* = 304 km 
s~^, ro — 0.16 kpc, /3 = 0.5, 5 = 10~^'^, which allows us to 
reproduce its surface brightness profiles both in the optical 
and in X-rays with a reasonable accuracy. In this case sinks 
and sources of mass and energy for the hot phase are locally 
compensated everywhere in the galaxy if the parameters a, 
b, and e* are constant. Hydrostatic equilibrium requires an 
adjustment of the temperature of the hot gas Tcq{b, To) and 



of the stellar velocity dispersion a, . This is achieved by ad- 
ditional restrictions applied to the values of b and To coupled 
by 



-Teq(b,To) =2a.^ 



(9) 



where TZ is the gas constant and th e mean molecular weight 
p — 0.62 [see (Kritsuk et al. 1998) for more deta ill. For a 



specific stellar rnass loss rate a = 5 x 1 s (Faber & 



Gallagher 197(: ; Sarazin fc White 1987 ) , hydrostatic equi 



librium is achieved with a reasonable mass deposition effi- 
ciency of b = 0.48 and a characteristic temperature of the 
heat source of To = 5 x 10^ K. This temperature is lower 
than the 7 x 10^ K corresponding to a SN Type la rate of 
O.55/175 SNu (van den Bergh & Tammann 1991), but it is 
higher than th e 3 x 10 K obtained fo r a supernova rate of 
0.18 ;i?5 SNu dCappellaro et al. 1999[ ). 

Note, that we do not claim NGC 4472 to be a perfect 
galaxy with an isothermal hot ISM, an isothermal stellar 
system and an isothermal dark halo. Rather we construct 
a simple model whose physical properties are based on ob- 
servations of NGC 4472. The use of a restricted set of con- 
trol parameters allows us to investigate in which way(s) our 
model differs from the real object. 



3 NUMERICAL METHOD 

The evolution ary equations arc solv ed by means of the 
AMRA code ( |Plewa fc Miiller 200c| ), which is based on 
the direct Eulerian version of the PPM method of Colella 



& Woodward (1984) on an adapti vely refined hierarchica l 
system of completely nested grids (Berger & Colella 1989). 



For the present ID and 2D axisymmetric simulations we 
use spherical coordinates. The computational domain ex- 
tends from Tin = 100 pc (rin — pc in 1-D simulations) to 
^out ~ 400 kpc in radius and covers (in 2D simulations) the 
wedge 10° < 9 < 170°. Periodic boundary conditions are 
imposed in the angular directionjT] In the radial direction 
hydrostatic equilibrium is assumed at rin and at rout, i.e. the 
radial (and theta) velocity is set to zero and the boundary 
pressure is defined in such a way that the pressure gradient 
at the boundary is equal to the hydrostatic one. We allow 
for up to 4 levels of refinement with density and pressure 
jumps of 10 per cent used as refinement criteria. With re- 
finement ratios of 4 and 2 in radius and angle, respectively, 
the resolution of our 2-D models is equivalent to that of a 
uniform grid run with 32768 x 160 zones (Ar ~ 12 pc and 
A6 = 1°). The source terms are calculated explicitly (Plewa 



1995 ) using the operator split approach. 



4 RESULTS 



t Unlike reflecting, periodic boundary conditions allow for a free 
gas flow across the sides of the computational volume, wliich is 
closer to the 3-D case. Switching from periodic to reflecting con- 
ditions in the steady convective regime does only modify the ex- 
treme values of the velocities. The main flow remains basically 
unchanged. 
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Figure 1. Quasi-steady state inflow (initial density perturbation e = 0.1). The left panel shows the density, temperature and velocity 
distributions at i = 0, 50, 100, 150, 200, 250, and 300 Myr for an inflow solution obtained with intermediate spatial resolution of 12 pc. 
The right panel shows hydrodynamic variables at t = 60 Myr. The solutions obtained with a spatial resolution of 3, 12, and 48 pc are 
shown by solid, dashed, and dotted lines, respectively. The heavy lines refer to the initial conditions. 



4.1 Spherically symmetric solutions 

In the presence of heat conduction the equilibrium rec ycling 



model for the h ot gas in giant ellipticals is meta-stable ( Krit- 
suk et <1. 1998). It means that, while arbitrary infinitesimal 



perturbations do not grow, perturbations of a finite ampli- 
tude can switch the system to a different equilibrium state 
or to a nonequilibrium catastrophic evolutionary regime. 

Since we are mostly interested in studying thermal ef- 
fects and since the characteristic thermal time scale is usu- 
ally longer than the dynamical time, it is quite natural to 
examine the response of the system to a class of perturba- 
tions, which do not disturb the hydrostatic equilibrium. In 
general, the system can be perturbed in a number of differ- 



ent ways. However, the approximation of spherical symme- 
try, which allows us to reduce the problem to 1-D, is rather 
restrictive. 

The simplest perturbation of this class is a global den- 
sity perturbation 5p — epaqir) which, at first sight, may 
seem to be unphysical because of its global nature. How- 
ever, due to stratification of the hot gas, the thermal time 
scale depends strongly on radius, i.e. for a time scale of a few 
10* yr the perturbation only affects the inner region because 
the rest of the hot gas has not enough time to react. Hence, 
for a sufficiently short time interval the global perturbation 
is equivalent to a more local one. 



© 2001 RAS, MNRAS 000, §-0 



Convectively Unstable Outflows 5 



4-1.1 Cooling inflows 

We have computed a series of spherically symmetric so- 
lutions with initial global density perturbation amplitudes 
e — 0.1 and 0.2, and compared the results obtained with dif- 
ferent spatial resolution in the central cooling region. Left 
and right panels of Fig. |l| illustrate our results for e = 0.1. 
The left panel shows distributions of density temperature 
and velocity at i = 0, 50, 100, 150, 200, 250, and 300 Myr 
for our intermediate resolution run with a radial spacing of 
12 pc. The heavy lines indicate the initial conditions. The 
quasi-steady regime is already established at t ~ 50 Myr 
after having passed through the linear growth regime which 
covers most of the earlier epochs. Then, for a very short time, 
a transient dynamic inflow occurs. Afterwards the flow re- 
mains quasi-steady. The gas temperature in the cooling core 
drops below the equilibrium value Toq by a factor of ~ 10. 
The density distribution shows a central cusp p oc r'^ with 
a power index f = —0.6. Hence, this model predicts a com- 
pact central peak in the X-ray surface brightness. This peak 
can be observed if the instrument has sufficient angular res- 
olution, and if finite optical depth efi^ects in the inflowing 
dense core will not render the approximation of an optically 
thin plasma invalid. Note, that heat conduction and/or mag- 
netic pressure support can modify the stability of the initial 
equilibrium model and even sup press the formation of the 
cooling core (Kritsuk et al. 1998). 

The right panel shows the flow variables at t = 60 Myr 
soon after a compact cooling infiow has been established 
in the core region of ~ 1 kpc in radius. The effects of the 
finite grid resolution are most evident in the velocity plot. 
Changes in the grid resolution effectively modify the position 
of the inner boundary thus determining the lowest velocity 
and temperature, and the highest density values attained on 
the grid. 

At the coarsest grid spacing of 48 pc the compact in- 
fiow region (100 ^ r ^ 500 pc) is covered by less than 10 
grid zones and hence remains unresolved. Steep gradients of 
dependent variables introduce interpolation errors, and er- 
rors due to explicit calculation of the source terms. These 
errors caused by the finite difference scheme are amplified by 
the physical thermal instability and, as a result, for e — 0.1 
the numerical solution displays high frequency noise in the 
temperature distribution for t J> 65 Myr. Using the same 
spatial resolution but doubling the size of the initial pertur- 
bation (e = 0.2) makes the infiow even more dynamic and 
the noise amplitude becomes so large that the numerical 
solution eventually switches from the infiow to the outfiow 
regime (see Fig. El). 

Grid resolution determines the amount of numerical 
dissipation of PPM. Although PPM does not utilize artifi- 
cial viscosity explicitely, its interpolation scheme introduces 
some amount of dissipation, which depends on grid reso- 
lution. When increasing the resolution by a factor of 16 to 
Ar = 3 pc the decreased numerical dissipation is found to be 
too small to prevent the development of small scale conden- 
sations due to the physical thermal instabili ty inh erent to 
fiows with radiative cooling [see, e.g.. Hunter ( |l970| ), Binney 
& Tabor (1995)]. As a result, a quasi-steady infiow solution 
is obtained (solid lines in the right panel of Fig. |l|), which be- 
comes unstable soon after t = 200 Myr causing the breakup 




o 



Figure 2. Sketch of the two ahernative sohitions for the meta- 
stable hydrostatic equiUbrium in spherical symmetry: a cooling 
infiow (left) and a subsonic outflow (right). The regions of intense 
mass deposition are shown in sohd black in the vicinity of the 
centre in case of inflow and a condensation front in case of outflow. 



of the central density enhancement into multiple spherical 
condensation fronts. 

At an intermediate spatial resolution of 12 pc (left pan- 
els and dashed lines in the right panels of Fig. ^) the cooling 
infiow region is well resolved by ~ 32 grid zones. The numer- 
ical solution remains stable for at least 500 Myr beyond the 
initial relaxation to a steady state infiow at ~ 55 Myr. Note 
that at higher spatial resolution (solid lines; right panels) 
the profiles are smoother than at lower grid resolution (left 
panels). The 'ripples' present at lower resolution can be elim- 
inated, for instance, by including a small phenomeno logica l 
diffusion term in the equations, as in Binney & Tabor (1995), 
or by tak ing into account thermal conductivity [Kritsuk et 
al. ( |l998| )]. 

The numerical experiments discussed below are all per- 
formed with a minimum radial grid spacing of 12 pc. At a 
much coarser grid important hydrodynamic scales remain 
unresolved, whereas a much finer grid is no longer consis- 
tent with our basic assumption that the sinks and sources 
of mass and energy are distributed continuously. 

Preliminary results obtained in two dimensions assum- 
ing axial symmetry indicate that cooling infiows are stable 
against random small velocity perturbations. The stability is 
maintained due to advection of linear perturbations towards 
the inner boundary by the accelerating gas infiow. 



4-1-2 Subsonic outflows 

As an alternative to quasi-steady cooling inflows we have 
studied a class of subsonic outflow solutions, which can be 
obtained through negative global perturbations added to the 
equilibrium density. In such a case heating and mass supply 
from evolved stars in the galaxy overtake cooling and mass 
dropout, and instead of a compact cooling core a hot bubble 
develops in the centre (Fig. |^ . 

Fig. ^ shows an example of a subsonic outfiow solution 
at t = 0, 100, 200, and 300 Myr for e = -0.1. In addition to 
the density, temperature, and velocity distributions we also 
show the entropy, which is defined here as 



s = Cv In 



-f const. 



(10) 



Numerical values for the entropy s axe given assuming 
const = 0. 

As the gas in the hot bubble expands it compresses the 
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ambient medium and thus creates a leading spherical quasi- 
isobaric condensation wave. The wave front is clearly seen 
as a peak in density and as a sharp minimum in the tem- 
perature profiles. The gas temperature in the bubble is a 
factor of ~ 2 higher than the equilibrium temperature Teq. 
It eventauUy saturates at a value of about 3 x 10^ K which 
is fixed by the heating rate at the centre. Typical expansion 
velocities in the bubble are about 10 km s~^. The maximum 
velocity continuously grows as the bubble expands while less 
and less dense outer layers of the hydrostatic corona start 
to feel the energy imbalance due to the initial global density 
perturbation. The bubble will expand beyond 300 Myr (the 
state shown in Fig. ^) because of the special form of the ini- 
tial perturbation. Thus, the solution will not tend towards a 
steady galactic wind. Since on a timescale of a few Gyrs the 
heating rate and stellar mass loss vary, such steady wind so- 
lutions are unrealistic, anyway. Negative velocities observed 
at r > 50 kpc are related to a sound wave, which is triggered 
at the beginning of the simulation by small deviations of the 
initial model from an exact hydrostatic state. 

The outflow solution is Rayleigh- Taylor unstable since 
the initially stable entropy stratification is modified by the 
expanding bubble. In particular, one may expect that the 
unstable local entropy inversion associated with the form- 
ing condensation front will be smeared out by convective 
motion, if the restriction of spherical symmetry of the fiow 
is relaxed, and, thus, such a sharp front will never form in 
reality. 



4.2 Unstable axisymmetric outflows 

In order to initialize the fiow in the axisymmetric case we 
impose besides a e = —0.1 global density perturbation also 
a small random velocity perturbation with an amplitude 
of 10^'' of the adiabatic sound speed. This additional ax- 
isymmetric perturbation helps us to initiate two-dimensional 
fiow. In the central region, 0.1 kpc < r < 4 kpc, the grid res- 
olution is 320 X 160 zones. The evolution has been followed 
for 300 Myr. 



4-2.1 Flow structure and evolution 

Convection starts to develop after t « 60 Myr from the be- 
ginning of the simulation in the form of small-scale vertical 
finger-like structures. The fingers originate the central region 
of unstable entropy inversion. High entropy material fiows 
upwards and low entropy material - downwards, restoring a 
marginally stable stratification with an adiabatic tempera- 
ture gradient 

?=fl-iU^. (11) 
dr \^ 7 y P dr ^ ' 

As the fingers flowing in opposite directions accelerate, 
Kelvin-Helmholtz roll-up transforms them into mushroom- 
l ike structures ty pical for the developing convective pattern 
( [nogamov 199S|, Section 6). 
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Figure 4. Four snapshots of the temperature field in the innermost 4 kpc region (from left to right) taken at t = 150, 200, 250, and 
300 Myr, respectively. The vertical scale along the symmetry axis shows the radial distance in kiloparsecs. Note that the left sectors in 
both panels [t = 150 and 250 Myr) are mirrored with respect to the right ones. The hot bubble at the centre can be seen in grades of 
red and the leading condensation front smeared by convection in blue. 




Figure 5. Same as Fig. 4 but for the entropy field. Gas motions triggered by the convective instability tend to level the inner part of 
the entropy distribution. A characteristic flow pattern in the core consists of ascending and descending mushroom-like structures. 

By t = 100 Myr a convective core with a radius of involved into the convective motion, too. As the temperature 

1 kpc forms and large scale structures become dominant. At in the hot bubble grows and the convective core expands, 

that time the radial velocities of the convective motion reach large scale convective structures with higher flow velocities 

±100 km s^^, and convection approaches a quasi-stationary appear, 
regime. Still later, when the outer gas layers start to feel 

the imbalance of mass and energy gains and losses due to At t = 300 Myr the hnear size of the largest mushroom 

the global initial density perturbation, these layers become is about 3 kpc and typical flow velocities are ^ 150 km s"\ 

A little bit earlier, at t = 280 Myr, the maximum Mach 
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number of 0.54 is recorded for a stream of upflowing hot gas. 
Thus, the convection driven by initially smalj^ (e = —0.1) 
deviations from local thermodynamic equilibrium remains 
subsonic. In a realistic 3-D situation one could expect even 
lower convective velocities due to the extra dimension avail- 
able for a redistribution of the kinetic energy. 

Snapshots of the temperature, entropy and vorticity dis- 
tributions for the inner 4 kpc are shown in Figs ^|-|^. The 
temperature maps clearly exhibit the expansion of the hot 
bubble and the distorted analogue of the spherically sym- 
metric leading condensation front, which is formed by the 
uppermost parts of the mushroom caps. The mushrooms 
are best seen in the entropy maps owing to the formation 
of an entropy 'plateau' in the core. The mean value of the 
entropy in the core steadily grows. Consequently, the con- 
vective pattern changes from blue via green and yellow to 
red. The vorticity maps highlight small scale structures in 
the velocity field and help in tracing the mushroom clouds 
as counter-rotating eddies forming the caps (Fig. The 
convective core as a whole has a well defined sharp outer 
boundary defined by the largest mushroom caps. The global 
maximum of the hot gas density is typically slightly off cen- 
tre. 



t Note, that a negative perturbation of 10 per cent in the gas 
density implies an about 20 per cent lower luminosity than in 
equilibrium, i.e. 20 per cent 'extra' heating. This corresponds to 
a net continuous supply of heat with a power of ~ 10'"^ ergs s~^. 
A conservative estimate of ~ 10*'^ ergs s~^ for the typical power 
supplied by AGNs in the form of mechanical energy of jets is mu ch 
larger than this value (see, e.g., Owen, Eilek & Kassim (2000)). 



For further details on the flow evolution we re- 
fer to a set of computer animations available at 



http : //www.mpa-garching.mpg. de/Hydro 



The typical power required to sustain radio core sources in el- 
liptical galaxies is about ~ 10'^^ ergs s~^. Our simulations show 
that a smaller (e = —0.01) initial density perturbation develops 
a convective core with a radius of 1.3 kpc by 300 Myr. 



4-2.2 Statistical properties: angular averages 

Analysis of the angular averages of the hydrodynamic vari- 
ables allows one to compare the properties of the new equi- 
librium!^ established in the convective core with those of the 
hydrostatic recycling model, and those of the outflow solu- 
tion. 

Fig. shows angular averages at 300 Myr for the central 
4 kpc region. Convection has smeared out the condensation 
front, which forms in the spherically symmetric case. As a 
result, the averaged convective solution falls in between the 
equilibrium model and the 1-D outflow. 

The gas density near the very centre is by a factor of ~ 3 
smaller than in the equilibrium, but still larger by a factor of 
~ 1.5 than in the 1-D case. Unlike in a compact inflow, there 
will be no sharp peak in the X-ray brightness in the central 
1 kpc in this case (see also the emissivity distribution in the 
upper panel of Fig. An observation with sub-arcsec res- 
olution of a nearby elliptical galaxy, like NGC 4472, would 
be a critical experiment to distinguish between the compact 
inflow and convective outflow realizations in the core. The 
relatively small density fluctuations about the mean indicate 
that compressibility is not very important, which is consis- 
tent with a subsonic convection regime. 



3 We refer to the term 'equilibrium' here in a statistical sense, 
implying that the moan convective flow is n u a si-steadv on time 
scales of a few 10 Myr. According to Fig. |l2| , the flow evolves 
slowly on a time scale of hundreds of Myrs which is longer than 
the ty pical life span of a convective element (See Sections 1.2.3| , 
4.2.5|). 
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Figure 7. Various hydrodynamic variables are shown as a function of radius for r < 4 kpc at 300 Myr. Heavy solid lines show angular 
averaged values, solid lines show max and min envelopes (i.e., the maximum/minimum angular value at the respective radius), dashed 
lines give the initial conditions, and dotted lines show the spherically symmetric solution at 300 Myr. The pressure is normalized to the 
initial hydrostatic value. The gas temperature is the emission measure weighted angular mean, while the velocity components (vertical 
u and horizontal v) and the entropy are mass weighted. 



The average temperature grows by a factor of ~ 2 in 
the central 1 kpc and has a shallow minimum at the po- 
sition of the condensation front. Local relative temperature 
variations in the convective flow are about 1.5, mainly at the 
interface between the convective core and the ambient hy- 
drostatic gas. The average pressure profile shows clear devia- 
tions from the hydrostatic stratification. Local pressure fiuc- 
tuations are < 5 per cent, and hence imply a quasi-isobaric 
regime of convection. The angular averaged entropy distri- 
bution is almost flat within a 2 kpc core and rises further 
outwards. There is no entropy inversion as in the 1-D case. 



i.e. the system has settled to a new equilibrium state via 
nonlinear saturation of the convective instability. 

Both vertical and horizontal gas velocities are quite 
high compared to typical cooling fiow values of about 
—20 km s^^. However, the mean radial velocity at t = 
300 Myr lies in the range —1.0 < u < 5.2kms~^, which is 
similar to the 1-D case, where Umax(< 4kpc) = 7.2 kms^^. 
Thus, there is no considerable mass flux in the radial direc- 
tion. 

In order to assess the efficiency of radiative cooling in 
the turbulent core, we define the volume emissivity e as 

e = Ptie + /5^A. (12) 
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The two terms on the rhs in the above equation describe 
radiative losses from the cooUng condensations and direct 
radiative cooUng of the hot gas phase, respectively. The up- 
per panel in Fig. |^ shows the distribution of the emissivity 
with radius for the inner 4 kpc. The peak central emissivity 
in the convective core is an order of magnitude lower than in 
hydrostatic equilibrium. For r > 1 kpc the emissivity distri- 
bution roughly follows the one of the recycling model. Local 
variations of the emissivity by a factor of ~ 3 are detected 
in the condensation front region between 2 and 3.5 kpc. 

The lower panel in Fig. ^ shows the relative rates of mass 
supply from stellar winds and mass deposition due to local 
thermal instabilities. Notice that for r < 1 kpc the mass sink 
is switched off due to a stabilization by intensive SN-heating 
(see equation |2|). Mass dropout dominates on average at radii 
from 2 to 3 kpc. The oscillations of the mass dropout rate are 
rather intense between 2 and 3.5 kpc, where the convective 
flow turns over. 

All of the best X-ray observations of elliptical galaxy X- 
ray halos show positive temperature g radients in the central 
>10 kpc ( |Brighenti fc Mathews 1997 ). XMM-TVetuton obser- 
vations of M87 and its X-ray halo report T ~ 1.3 keV in the 
central circle of 0.25' (~ 1.2 kpc) in radius, a 'steep' tem- 
perature rise to ~ 2.3 keV at ~ 3' (~ 15 kpc), and a shallow 
rise to ~ 2.9 keV in the range from 3' to 13' (Bohringer et 
al. 20o||). Yet the measured temperature gradients do not 
give us much information about the temperature distribu- 
tion within the central kpc. Only the latter would help us to 
verify the model assumptions concerning the thermal bal- 
ance of the hot gas near the centre. The measured gradients 
can rather point to a slight non-isothermality of the underly- 
ing gravitational potential, which traces the transition from 
the dark matter hal o /cluster dominated oute r region to th e 
stellar core of M87 (|Nulsen fc Bohringer 1995[ [Nulsen 1998|). 



At the same time, the flat (on the average) entropy dis- 
tribution in the core region of our model is supported by the 
assumed spherical symmetry of the underlying gravitational 
potential, average gas density distribution, and the net heat 
source. As we mention in the introduction, this assumption 
is restrictive. Even a small asymmetry, e.g., in the central 
heating efficiency would create a hybrid in/out-flow with a 
more stable (on the average) entropy stratification, i.e. with 
a mild positive average temperature gradient. As soon as the 
asymmetric heating is sufBcient and, still, centrally concen- 
trated (as in the case of AGN heating), such a flow would 
not necessarily display a central density/X-ray brightness 
cusp. 

Observationally, one can expect to see convection in the 
form of uprising hot bubbles interacting with a background 
flow of cooling gas. It is possible, that flows of this kind 
are observed in X-rays by Chan dra and in radi o on VLA in 



NGC 1275 ( [Fabian et al. 2000D and in M87 ( pwen et al 
200C|)7^uch flows seem to be largely subsonic, and the in- 
terface of a hot bubble with the background flow does not 
imply shoc k wav es unlike in the model of Heinz, Reynolds & 
Begelman (1998). Rather, it would be physically similar to 



the interface of our convective core with the surrounding hy- 
drostatic corona, displaying lower temperatures, larger mass 
deposition rates and a higher X-ray surface brightness. We 
shall consider such flows in detail elsewhere. 




12 3 4 

r (kpc) 

Figure 8. Mass and energy supply rates for the region r < 4 kpc 
at 300 Myr. The upper panel shows the volume emissivity. The 
lower panel displays the relative rates of mass supply to and 
dropout from the hot phase. Heavy solid lines show the angular 
averaged mean values, solid lines are the max and min envelopes 
(i.e., the maximum/minimum angular value at the respective ra- 
dius), and dashed lines give the initial conditions. 



4-2.3 Statistical properties: power spectra 

A spectral decomposition of a flow helps to reveal informa- 
tion which is not apparent otherwise, like the distribution of 
energy over length scales and effects of the finite grid resolu- 
tion. In order to reveal this information we have computed 
one-dimensional Fourier transforms along the ^-coordinate 
direction at various radii. Strong gravitational stratification 
of the gas density and the lack of periodicity in the radial di- 
rection would dominate the spectrum in the radial direction, 
and thus would complicate the analysis. 

For a scalar quantity of interest, q, we define the hori- 
zontal Fourier transform to be 



Q{r,k,t) 



1 



q{r,e,t)e 



-27rifcfl/AS 



(13) 



which is a function of time, radius, and the horizontal 
wavenumber k. Here AO is the horizontal periodic length. 
The power spectrum is the square of the modulus of its 
complex Fourier transform, namely | Q{r,k,t) f. We have 
used simple periodo grams [see, e.g.. Press, Teukolsky, Vet- 
terling & Flannery (1992)] as the power spectra estimates 
for discrete data sets given on the finest level grid covering 
the inner convective region of the computational domain. 
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Figure 9. Time averaged power spectra estimates of the vertical (-u, left panel) and horizontal (n, right panel) component of the gas 
velocity. The power spectra are shown for radii of 0.71 (solid line), 1.32 (long dashes), 1.93 (dash dot), 2.54 (dashed), and 3.15 kpc 
(dotted), respectively. The power spectra are averaged over 80 consecutive states covering the time interval 275-300 Myr. The spectra 
show no significant radial dependence. 



For the spectral analysis we have used the velocity field. 
When the density is constant, the power spectrum of the 
velocity field is similar to a kinetic energy spectrum. As 
the gas density varies only slightly about its mean value 
with angle at a given radius (see max and min envelopes 
in Fig. the density factor is always negligible. However, 
since radial modes are ignored in our ID analysis, the sum of 
the power spectra of the velocity components is not expected 
to be a measure of the total kinetic energy. 

Since instantaneous power spectra are noisy and since 
they fluctuate considerably in time, we use temporal av- 
erages to estimate the typical power distribution with 
wavenumber. The convective core expands as the source 
terms associated with the stellar system drive the hot bubble 
at the centre. Hence, we had to take averages over a period 
of time shorter than the expansion time scale (~ 10* yr) but 
longer than the lifetime of the largest mushroom-like struc- 
tures in our convective flow (~ 20 Myr, see Section 4.2.5). 
Averaging over a time interval 25 Myr sufficiently suppresses 
the oscillations in the mean spectra at small wavenumbers. 
We use angle brackets < . . . > in the following to denote the 
computed temporal averages. 

Fig. 1^ shows the power spectra of the two velocity com- 
ponents at radii of 0.71, 1.32, 1.93, 2.54, and 3.15 kpc, re- 
spectively. Obviously, the power is concentra ted to wards the 
large scales, as already mentioned in section 4.2.1 when dis- 



cussing the settling of the convective regime. The spectra do 
not show any significant radial dependence. The mean slope 
p of the spectra <| Q p>oc is p ~ — |. Following the 
work of Porter & Woodward (1994), who investigated the 



numerical dissipation of the PPM scheme, we estimate that 
in our simulations numerical dissipation starts to become 
important for k > 10, i.e. for structures smaller than about 
16 angular zones. 

Batchelor (1969) in his paper on two-dimensional 
isotropic turbulence derived an analogue of the Kolmogorov 
energy spectrum. For 2D flows the spectrum of vorticity 
I < (V X v)'^ > is universal. If Q{k) is the vorticity spec- 
trum, then Q{k) oc is the analogue of Kolmogorov's 



power law. For the energy spectrum of 2D flows one has 
E{k) cx k''^Q{k) cx fc"^. From our Fourier analysis we ob- 
tain as expected |[/'^| cx k""^ and [V'^l oc k~^ in the interval 
k £ [3,30]. This explains the transition from small scale to 
large scale mushrooms when convection sets in. 



4.2.4 Average transport properties 

The characteristics of the convective flow in the inner region 
can also be examined by calculating the angular averaged 
energy and a set of energy fluxes, which we d efine below 
following Hurlburt, Toomre & Massaguer ( |l986| ). 

The angular integrated total energy equation (||) can be 
expressed as 



S + + Fc + Fk) = S, 

at r^ or 



with 

E = 2tt 



ErsmOde, 



(14) 



(15) 



where F^ is the heat flux due to a spherically symmetric 
expansion or contraction of the model, Fc is the 'enthalpy' 
or convective heat flux, _Fk is the kinetic heat flux, and S 
is the net energy source term. These fluxes and the source 
term are deflned as 



Fc = 27r 



271 



up ■ 



up ■ 



up ■ 



- + — 7- \ r smO dd, 

p 2 dr ' 



e + 



— r^'smOde, 



ap^e, — pti 7e + 



(16) 



(17) 



(18) 



'sin 6'd6l.(19) 
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Figure 10. Total angular integrated energy flux _Ft (upper panel) 
and net energy source S (lower panel) at t = 10, 30, 50, and 
70 Myr (lines from top to bottom in each panel), respectively. 



The above definitions involve spatially fluctuating variables 
/ = / + /', where / is the angular mean value and /' is the 
fluctuation about that mean value. They allow us to dis- 
cuss correlations between the radial velocity and a number 
of fluctuating quantities. Such correlations provide useful 
insights into energy generation and energy transport mech- 
anisms in the convective interior of the hot ISM. 

We further define a total energy fiux Ft = F^+Fc+F^. 
In a steady state with S = Q the total flux does not de- 
pend on radius r. If a convective flow achieves a statistically 
steady state, then for sufficiently long time averages the av- 
erage total energy fiux 



<Ft> = <Fe> + <Fc> + <Fk> 



(20) 



does not dependent on radius as well, while the individual 
contributing fiuxes may vary considerably with r. Before 
turning to the fluctuations let us consider the inner 4 kpc re- 
gion where convective motion occurs in a mean 'background' 
flow after the onset of the instability. Fig.|l^ shows the ra- 
dial dependence of the total energy flux Ft and net source 
S at four representative moments of time before and shortly 
after the instability sets in. At these times the total fiux is 
directed outwards and increases from zero at r = to values 
around 2.5 x 10^" erg/s at r = 4 kpc mainly due to a slow 
overall expansion of the gas. This expansion is stimulated 
partly by an outgoing weak sound wave generated by the 
initially imperfect hydrostatic equilibrium, and also by the 
increasing source terms blowing a hot bubble at the centre. 
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Figure 11. Mean energy fluxes < F > (upper panel) and mean 
net energy source < S > (lower panel) averaged over 295.5 < 
t < 308.3 Myr. The heavy solid line shows the mean total flux 
< -Ft >, and the solid, dotted, and dashed lines show the heat 
flux due to a spherically symmetric expansion or contraction of 
the model < _Fe >, the convective flux < Fq >, and the kinetic 
flux < _F]j >, respectively. 



Since the fiow is spherically symmetric for t < 50 Myr, the 
major contribution to the total fiux is due to gas expansion 
in an external gravity field, i.e. Ft — Fe- The small wrinkles 
visible at r ~ 1 kpc in the flux diagram at t = 70 Myr hint 
towards a contribution from fiuctuations of the developing 
convective fiow, but still Ft ~ Fe. By t = 300 Myr, when 
the convective fiow already covers the whole inner region 
(r < 3 kpc), the global expansion driven by the action of 
the source terms remains the dominant contributor to the 
total energy fiux (Fig. ^l|). Note that the mean energy fiux 
carried away from the inner region aX t ^ 300 Myr amounts 
to about 7 per cent of the luminosity of the hot gas in the 
whole computational domain r < 400 kpc (see Fig. [l^ ). 

Our analysis indicates that the fiuctuations play only 
a secondary role for the energy transport. The mean angu- 
lar integrated convective fiux, < Fc >, is positive every- 
where except for the region around r ~ 2.3 kpc, where the 
mean source < S > attains its minimum, yielding an up- 
ward convective transport of thermal energy. Like in case of 
an incompressible fiuid a positive fiux value corresponds to 
motions where hot fluid parcels rise and cold ones sink. The 
depression in < Fc > at r « 2.3 kpc spans the same range in 
radius as the condensation front smeared by convection. In 
this region mass dropout due to thermal instabilities domi- 
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Figure 12. Temporal behaviour of the luminosity L (thin curve) 
and net mass deposition rate M (thick curve) integrated over the 
whole computational domain. 



nates over mass supply from stellar winds, and the enthalpy 
flux, carried by fluctuations from the inner region, decreases. 

The mean kinetic flux < Fk > steadily grows in abso- 
lute value, as the convective core becomes more extended 
giving rise to larger convective structures with higher ve- 
locities. The kinetic flux is positive (i.e. directed outwards) 
at r < 2.5 kpc and weakly negative further out in the re- 
gion where convective motions overshoot through the con- 
densation front. The sign of Fk is related to differences in 
strength and filling factor between upward and downward di- 
rected motions. In our two-dimensional simulations upflows 
are considerably faster than downflows in the interior of the 
convective core and there is no apparent asymmetry in the 
filling factors (see the angular velocity diagram in Fig. ^). 



4-2.5 Integral characteristics 

In order to follow the evolution of global characteristics of 
the flow we computed the total hot gas luminosity L and the 
net mass deposition rate M as integrals over the whole com- 
putational domain (Fig.^). The computed total luminosity 
is close to the X-ray luminosity Lx ~ 7.2 x 10*^ erg s~^ of 
NGC 4472 at a distance of 17 Mpc. There is a mild tendency 
for the luminosity to grow with time as the gas shed by stars 
in the galaxy slowly accumulates in the convective core. Low 
amplitude oscillations in the luminosity on a time scale of 
r ~ 20 — 30 Myr correspond to a life cycle of the largest 
convective structures. Assuming a linear scale of I ~ 3 kpc 
and r — 25 Myr, one can estimate a characteristic velocity 
1/t ~ 120kms~^. This value is typical for the convective 
core which developed in the simulations by t ~ 300 Myr (see 
Fig. ^) . A time scale of 25 Myr is also a reasonable choice 
for the time averaging interval to be used when computing 
mean power spectra. 

The net mass deposition rate M, given by the volume 
integral 



M 



(apt - pti)dV , 



(21) 



is initially positive, since a negative global density perturba- 



tion Sp/p = —0.1 breaks the balance between the mass sup- 
ply to the hot phase from stellar winds and mass loss from 
the hot phase due to local thermal instabilities. The mass 
supply from stars always dominates, and M grows from 0.5 
to 0.75 Meyr-^ (Fig.||) since the system evolves further 
away from hydrostatic equilibrium. For our model galaxy 
with a stellar mass M« « 8.5 x 10^^ M0 the rate of steady 
supply of gas from stars, a,Af» ~ 1.2 Mq yr~^, implies that 
the gas condensation rate varies from 0.7 at f = to 0.45 Mq 
yr^^ at 300 Myr. As one can see in Fig. |l2[ M slightly os- 
cillates in phase with the X-ray luminosity. Maxima of M 
record moments when streams of rising hot material form 
new large scale mushroom-like structures. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have studied hot gas flows originating in 
the central few kpc of hot galactic coronae in response to a 
slight energy supply imbalance in a model initially in equi- 
librium. It is shown that a negative energy budget drives 
a quasi-steady compact cooling inflow, which is stable to 
small perturbations. Excessive heating creates an isentropic 
(on average) convective core in the vicinity of the galaxy cen- 
tre. Provided the net energy gain in the core is not too high, 
the convection remains subsonic with typical flow veloci- 
ties (both, radial and angular) of the order of ±150 km s^'^. 
The characteristic convective pattern consists of a manifold 
of mushroom-like structures of different scales. The power 
spectra of the velocity components are dominated by the 
largest-scale mushrooms. The emission weighted gas temper- 
ature in the core is slightly rising towards the centre. Unlike 
in the case of a cooling inflow, the X-ray surface brightness of 
a convective core does not display a sharp maximum at the 
centre. The inner parts of the convective core are dominated 
by mass supply from evolved stars, while in the colder outer 
shell mass drop out dominates the mass budget. Large-scale 
mushrooms (mainly their caps) can probably be marginally 
detected by high resolution X-ray images with good photon 
statistics, as their volume emissivity is about a factor of 2 
higher than that of the surrounding gas. 

We have considered an idealized case with a spherically 
symmetric heat source and small deviations from 'thermal' 
equilibrium near the centre. However, we can extrapolate 
our results to a more realistic hydrodynamic picture. First, 
stronger heating in the vicinity of the centre can drive a su- 
personic convection with velocities in excess of 400kms~^ 
and produce a much larger (here 3-4 kpc) convective region. 
Second, a Type la SN heating rate which is lower all over the 
galaxy, and an excessive heating by an AGN in the central 
region will create a hybrid flow with an inner convective core 
and an outer 'cooling flow', which will vanish in convection 
when the inflow velocity is very subsonic. Our preliminary 
results indicate that such a hybrid dynamic flow can be sta- 
ble and can last as long as the central source is active. Third, 
a low Type la SN heating rate without extra energy supply 
near the centre will drive a more extended massive cooling 
inflow with an even more pronounced central surface bright- 
ness peak. It would be harder to turn such a dynamic inflow 
into an outflow by a strong perturbation (e.g., a burst of 
AGN activity, or a merger) than a weaker compact inflow. 
Finally, we expect that an asymmetric central heat source 
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in excess of the equilibrium value would create bubble(s) 
filled with a (mean) isentropic convective flow with proper- 
ties similar to the ones described above. 

Chandra, Newton and future X-ray missions will teach 
us more about the hydrodynamics of 'cooling flows'. 
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